Effects of color superconductivity on the structure and formation of compact stars 



O 

o 

(N 
C 

a 



> 

in 
o 

o 

o 



o 

=5 



X 



A. Drago a , A. Lavagno 6 , G. Pagliara" 
1 Dipartimento di Fisica, Universita di Ferrara and INFN, Sezione di Ferrara, 44100 Ferrara, Italy and 
b Dipartimento di Fisica, Politecnico di Torino and INFN, Sezione di Torino, 10129 Torino, Italy 

We show that if color superconducting quark matter forms in hybrid or quark stars it is possible 
to satisfy most of recent observational boundaries on masses and radii of compact stellar objects. 
An energy of the order of 10 53 erg is released in the conversion from a (metastable) hadronic star 
into a (stable) hybrid or quark star in presence of a color superconducting phase. If the conversion 
occurs immediately after the deleptonization of the proto-neutron star, the released energy can 
help Supernovae to explode. If the conversion is delayed the energy released can power a Gamma 
Ray Burst. A delay between the Supernova and the subsequent Gamma Ray Burst is possible, in 
agreement with the delay proposed in recent analysis of astrophysical data. 

PACS numbers: 25.75.Nq, 26.60.+C, 26.50.+X 



The new accumulating data from X-ray satellites pro- 
vide important information on the structure and forma- 
tion of compact stellar objects. Concerning the struc- 
ture, these data, are at first sight difficult to interpret 
in a unique and self-consistent theoretical scenario, since 
some of the observations are indicating rather small radii 
and other observations are indicating large values for the 
mass of the star. Concerning the formation scenario, 
crucial information are provided by the very recent ob- 
servations of Gamma- Ray Bursts (GRB), indicating the 
possibility that some of the GRBs are associated with a 
Supernova (SN) explosion, ft has not yet been clarified 
if the two explosions are always simultaneous or if, at 
least in a few cases, a time delay can exist, with the SN 
preceding the GRB QHEIHE). 

The effect of the transition to deconfined Quark Matter 
(QM) on explosive processes like SNs and GRBs has been 
discussed by many authors. In particular, the possibility 
that deconfinement takes place during the core-collapse 
of massive stars at the moment of the bounce, has been 
discussed e.g. in Refs.0, and this mechanism could 
help the SN to explode by increasing the mechanical en- 
ergy associated with the bounce. However, it seems more 
plausible that deconfinement takes place only when the 
proto-neutron star (PNS) has deleptonized and cooled 
down to a temperature of a few MeV 0, . The energy 
released in the conversion to QM produces a refreshed 
neutrino flux which can help the supernova to explode 
in a neutrino-driven scheme. Finally, another scenario is 
possible in which neutron stars having a small enough 
mass can exist as metastable Hadronic Star (HS) if a 
non-vanishing surface tension is present at the interface 
between Hadronic Matter (HM) and QM. The process of 
quark deconfinement can then be a powerful source for 
GRBs and it can also explain the possible delay between 
a SN explosion and the subsequent GRB [To| . 

In recent years, many theoretical works have investi- 
gated the possible formation of a diquark condensate in 
QM, at densities reachable in the core of a compact star 
[llL Il3 | . The formation of this condensate can deeply 
modify the structure of the star [3 13 13 EI- We 
present here an extension of the previous works, showing 



that it is possible to satisfy the existing boundaries on 
mass and radius of a compact stellar object if a diquark 
condensate forms in a Hybrid Star (HyS) or a Quark Star 
(QS). Moreover, we show that the formation of diquark 
condensate can significantly increase the energy released 
in the conversion from a purely HS into a more stable 
star containing deconfined QM. 

To describe the high density Equation of State (EOS) 
of matter we adopt standard models in the various den- 
sity ranges. Concerning the hadronic phase we use rel- 
ativistic non-linear models 

[3d At 

very low density 

we use the standard EOSs of Refs. [3 E3 . For the QM 
phase we adopt a MIT-bag like model in which the for- 
mation of a diquark condensate is taken into account. 
To connect the two phases of our EOS, we impose Gibbs 
equilibrium conditions. 

It is widely accepted that the Color-Flavor Locking 
phase (CFL) is the real ground state of QCD at asymp- 
toticly large densities. We are interested in the bulk prop- 
erties of a compact star and we adopt the simple scheme 
proposed in Refs. [313 where the thermodynamic po- 
tential is given by the sum of two contributions. The first 
term corresponds to a "fictional" state of unpaired QM in 
which all quarks have a common Fermi momentum cho- 
sen to minimize the thermodynamic potential. The other 
term is the binding energy A of the diquark condensate 
expanded up to order (A/^) 2 . In Rei.[l4j the gap is as- 
sumed to be independent on the chemical potential \i. In 
the present calculation we consider a /i dependent gap 
resulting from the solution of the gap equation • The 
resulting QM EOS reads: 
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with v = 2/i— y M 2 + ~~3~ > an d the quark density p is 
calculated numerically by deriving the thermodynamic 
potential respect to fj,. Pressure and energy density read: 
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FIG. 1: Mass-radius plane with observational limits and repre- 
sentative theoretical curves: thick solid line indicates CFL quark 
stars, thick dot-dashed line CFL hybrid stars, thick-dashed line 
hadronic stars (see text). Observational limits from: (a) Sanwal 
et al. 2002 0, (b) Cottam et al. 2002 0, (c) Quaintrell et 
al. 2003 IH, (d) Heinke et al. 2003 (e),(g) Dey et al. 1998 
[Sj, (f) Li et al. 1999 (h) Burwitz et al. 2002 @. 

E/V = n CF L{n)+W + B + n e {ii e ) + ii e p e . (3) 

In Fig.l we have collected most of the analysis of data 
from X-ray satellites, concerning: masses and radii of com- 
pact stellar objects j^HH |H \M IH Observing 
Fig.l, we notice that the constraints coming from a few 
data sets (labeled "e" , "f" j^jj "g" and maybe also con- 
straint "h" |40j ) indicate rather unambiguously the ex- 
istence of very compact stellar objects, having a radius 
smaller than ~ 10 km. At the contrary, at least in one 
case ("a" in the figure), the analysis of the data suggests 
the existence of stellar objects having radii of the order 
of 12 km or larger, if their mass is of the order of 1.4 M Q . 
We recall that it is difficult from an astrophysical view- 
point to generate compact stellar objects having a mass 
smaller than 1 M©. Therefore the most likely interpre- 
tation of constraint "a" is that the corresponding stellar 
object does not belong to the same class of objects which 
have a radius smaller than ~ 10 km. Concerning con- 
straint "b" , it can be satisfied both with a very compact 
star or with a star having a larger radius. The appar- 
ent contradiction between the constraints "e" , "f" , "g" 
and the constraint "a" can be easily accommodated in 
our scheme, since it can be the signal of the existence 
of mctastable purely HS which can collapse into a stable 
configuration when deconfined QM forms inside the star. 

Finally, constraints ("c"H3 and "d"H2 ) do not pro- 
vide stringent limits on the radius of the star, but they 
put strong constraints on the lower value of its mass. It 
is in general not easy to obtain stellar configurations hav- 
ing both large masses and very small radii. As we will 
see, the existence of an energy gap associated with the di- 
quark condensate helps in circumventing this difficulty, 
since the effect of the gap is to increase the maximum 
mass of stars having a huge content of pure QM. 

In Fig.l we show a few theoretical M-R relations which 
correspond to the scenario we are proposing. More pre- 
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FIG. 2: Gap as function of the chemical potential, for four dif- 
ferent parameter sets. 



cisely, we show a thick-dashed line corresponding to HSs 
(GM1), a thick dot-dashed line corresponding to HySs 
(GM1, B 1 / 4 = 170 MeV, A 2 ) and a thick solid line corre- 
sponding to QSs (B 1 / 4 = 170 MeV, AA Similar shapes 
can be obtained using the EOS of Ref.[l9|. Both the HyS 
and the QS lines can satisfy essentially all the constraints 
derived from observations. The shapes of the gaps Aj 
are shown in Fig. 2. In conclusion, in our scheme most 
of the compact stars are either HySs or QSs having a 
mass in the range 1.2 — 1.8M Q and a radius ~ 8.5 — 10 
km. Stars having a significantly larger radius (as the one 
suggested by constraint "a" ) correspond in our scheme to 
metastable HSs which can exist if their mass is not too 
large, as we show in the following. 

Let us now discuss AE, the energy released in the con- 
version from HS to HyS or QS. AE is the difference be- 
tween the gravitational mass of the HS and that of the 
final HyS or QS having the same baryonic mass. As 
mentioned in the introduction, a possibility is that de- 
confinement takes place a few seconds after the bounce, 
when the PNS has deleptonized and its temperature has 
dropped down 0, 0- In particular, for stars having a 
small mass the formation of QM takes place only at T < 
few MeV. Notice that for a star having a mass of order 
1.4M Q and using the relativistic EOSs discussed in this 
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TABLE I: Energy released AE (measured in foe=10 5 erg) in 
the conversion from a 1.4 Mq hadronic star into the hybrid or 
quark star having the same baryonic mass (labeled with a •), 
for various sets of model parameters. BH indicates that the 
hadronic star collapses to a Black Hole. A dash (-) indicates 
situations in which the Gibbs construction does not provide 
a mechanically stable EOS. 
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paper, hyperons are present in the initial configuration, 
since the typical mass at which hyperons start forming 
is ~ IMq. The energy released during quark deconfinc- 
ment powers a new neutrino flux which can be useful to 
make the supernova explode. AE is shown in Table I, 
for a PNS having a mass of 1.4 M Q . As it can be seen, 
AE can be as large as 10 53 erg, if the final configura- 
tion corresponds to an HyS and three times as large if a 
QS is obtained. The effect of the gap is to increase the 
energy released and to allow QS configurations in cases 
where an HyS would be obtained in the absence of quark 
pairing. Let us now remark that the deconfincment tran- 
sition can be delayed if a non-vanishing surface tension 
at the interface between HM and QM exists and if the 
mass of the HS is not too large. This possibility was 
not discussed in Ref.|9( and it is the main ingredient of 
our model. To compute the time needed to form QM 
we use the technique of quantum tunneling nucleation. 
We can assume that the temperature has no effect in our 
scheme because, as discussed above, when QM forms the 
temperature is so low that only quantum tunneling is a 
practicable mechanism. 

In Ref. [l0( it was proposed that the central density of a 
pure HS (containing hyperons) can increase, due to spin 
down or mass accretion, until its value approaches the 
deconfinement critical density. At this point a spherical 
virtual drop of QM can form. The potential energy for 
fluctuations of the drop radius R has the form [3l| : 



U(R) = -TrR 3 n q (fi q 



AnaR 2 + 8n~fR (4) 



where n q is the quark baryon density, /ih and \i q are the 
hadronic and quark chemical potentials, all computed at 
a fixed pressure P, and a is the surface tension for the 
interface separating quarks from hadrons. Finally, the 
term containing 7 is the so called curvature energy. For 
cr we use standard values from 10 to 40 MeV/fm 2 and we 
assume that it takes into account, in a effective way, also 
the curvature energy. The value of a was estimated in 
Ref. to be ~ 10 MeV/fm 2 . Values for a larger than 
~ 30 MeV/fm 2 are probably not useful at the light of the 
results of Refs. mil. 

The calculation proceeds in the usual way: after the 
computation (in WKB approximation) of the ground 
state energy Eq and of the oscillation frequency vq of the 
virtual QM drop in the potential well U(R), it is possi- 
ble to calculate in a relativistic frame the probability of 
tunneling as |35| 



p = exp[-2 



A(E ) 



(5) 



where 



R4 



A(E) = / dRy/[2M(R) + E- U(R)][U(R) - E] 

(6) 

Here M{R) — Anph (l — R 3 , Ph is the hadronic en- 
ergy density and rih, n q are the baryonic densities at a 



same and given pressure in the hadronic and quark phase, 
respectively. Finally, R± are the classical turning points. 
The nucleation time is then equal to r = {vqpqNc)^ 1 , 
where N c is the number of centers of dro plet formation 
in the star, and it is of the order of 10 48 35] . r can be 
extremely long if the mass of the metastable star is small 
enough but, via mass accretion, it can be reduced from 
values of the order of the age of the universe down to 
a value of the order of days or years. We can therefore 
determine the critical mass M cr of the metastable HS for 
which the nucleation time corresponds to a fixed small 
value (1 year in Tab. 1). 

In Table II we show the value of M cr for various sets 
of model parameters. In the conversion process from a 
metastable HS into an HyS or a QS a huge amount of 
energy AE is released. We see in Table II that the forma- 
tion of a CFL phase allows to obtain values for AE which 
are one order of magnitude larger than the correspond- 
ing AE of the unpaired QM case (A = 0). Moreover, we 
can observe that AE depends both on magnitude and 
position of the gap. 

In the model we are presenting, the GRB is due to the 
cooling of the justly formed HyS or QS via neutrino - anti- 
neutrino emission. The subsequent neutrino-antincutrino 
annihilation generates the GRB. In our scenario the du- 
ration of the prompt emission of the GRB is therefore 
regulated by two mechanisms: 1) the time needed for the 
conversion of the HS into a HyS or QS, once a critical- 
size droplet is formed and 2) the cooling time of the justly 
formed HyS or QS. Concerning the time needed for the 
conversion into QM of at least a fraction of the star, the 
seminal work by |36| has been reconsidered by |37l ] , where 
it has been shown that the stellar conversion is a very fast 
process, having a duration much shorter than Is. On the 
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TABLE II: Energy released AE in the conversion to hybrid 
or quark star, for various sets of model parameters, assuming 
the hadronic star mean life-time r = 1 yr (see text). M cr 
is the gravitational mass of the hadronic star at which the 
transition takes place, for fixed values of the surface tension 
a and of the mean life-time r. Notations as in Tab. 1 
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other hand, the neutrino trapping time, which provides 
the cooling time of a compact object, is of the order of a 
few ten seconds jHslj, and it gives the typical duration of 
the GRB in our model. 

In conclusion, comparing the theoretical mass-radius 
curves with recent observational data, we find that color 
superconductivity is a crucial ingredient in order to sat- 
isfy all the constraints coming from observations. The 
difficult problem posed by astrophysical data indicating 
the existence of stars which are both very compact and 
rather massive can be solved either with hybrid or quark 
stars. Concerning hybrid stars, the gap increases the 
maximum mass of the stable configuration, while keep- 
ing the corresponding radius < 10 km. 

The superconducting gap affects also deeply the energy 
released in the conversion from hadronic star into hybrid 
or quark star. We assume that the deconfinement transi- 



tion only takes place when the star has deleptonized and 
cooled down, in agreement with the results of Ref. 0,0. 
If deconfinement occurs immediately after deleptoniza- 
tion, the energy released can help the SN to explode. 
If, at variance, the transition is delayed, a metastable 
hadronic star can form. Its subsequent transition to a 
stable configuration, containing deconfined quark mat- 
ter, can power a GRB via the annihilation of neutrinos 
and antineutrinos emitted during the cooling of the newly 
formed compact star. The energy released is significantly 
increased by the effect of the chemical-potential depen- 
dent superconducting gap and it can reach a value of 
the order of 10 53 erg. The proposed mechanism could 
explain recent observations indicating a possible delay 
between a Supernova and the subsequent Gamma Ray 
Burst dSQ. 
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